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Abstract. Black hole spins affect the efficiency of the "classical" accretion 
processes, hence the radiative output from quasars. Spins also determine how 
much energy is extractable from the hole itself. Recently it became clear that 
massive black hole spins also affect the retention of black holes in galaxy, because 
of the impulsive "gravitational recoil", up to thousands km/s, due to anisotropic 
emission of gravitational waves at merger. I discuss here the evolution of massive 
black hole spins along the cosmic history, due to the combination of mergers and 
accretion events. I describe recent simulations of accreting black holes in merger 
remnants, and discuss the implication for the spins of black holes in quasars. 



1. Introduction 

Black holes, as physical entities, span the full range of masses, from tiny black 
holes predicted by string theory, to monsters weighing by themselves almost as 
much as a dwarf galaxy (MBHs). Notwithstanding the several orders of mag- 
nitude difference between the smallest and the largest black hole known, all of 
them can be described by only three parameters: mass, spin and charge. Astro- 
physical black holes are even simpler, as charge can be neglected. So, besides 
their masses, Mbhj astrophysical black holes are completely characterized by 
their dimensionless spin parameter, a = J hi J max = c Jh/G Mg H , where is 
the angular momentum of the black hole, and < a < 1. 

1.1. Radiative efficiency 

In radiatively efficient, geometrically thin accretion disks the mass-to-energy 
conversion efficiency, e, equals e = 1 — E/Mbhc 2 , where E is the binding energy 
per unit mass of a particle in the innermost stable circular orbit (ISCO). The lo- 
cation of the ISCO depends solely on a black hole spin, shrinking by a factor of 6 
between a non-rotating hole and its maximally rotating counterpart- The closer 
the ISCO is to the horizon, the higher the mass-to-energy conversion efficiency, 
which increases from 6% to 42% in the above example. The mass-to-energy 



1 We will use th e term "maximally rotating" for a black hole with a = 1, although we note that 
iThorne! (1 19741 ) showed that accretion driven spin-up is limited to a = 0.998. Magnetic fields 
connecting material in the disk and the plunging region may further reduce the equilibrium spin. 
Magnetohydrodynamic simu lations for a series of thick accretion disks suggest an asymptotic 
equilibrium spin at a m 0.9 ([Gammic ct al. 2004). The location of the ISCO depends also on 
the particle being on a prograde or retrograde orbit. The ISCO for a particle on a retrograde 
orbit around a hole with a = 1 is 9 times larger than for its prograde counterpart. 
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conversion directly affects the mass-growth rate of black holes: high efficiency 
implies slow growth. More precisely, for a hole accreting at the Eddington rate, 
the black hole mass increases with time as: 

M(t)=M(0)e*p(i^^-), (1) 

where ^Edd = 0.45 Gyr. The higher the spin, the higher e, implying longer 
timescales to grow the MBH mass by the same number of e-foldings. Going from 
e = 0.06 to e = 0.42, the difference in M(t) amounts to 6 orders of magnitude, 
at t = tEdd- The typical spin therefore affects the overall mass-growth of MBHs, 
and the duty cycle of quasars. 

The r adiative effici ency is also the fundamental free paramete r for the Soltan 



argum ent (jSoltanll 19821 ) and, more recently, synthesis models (e.g. . Merloni fc Heinz 



2008) which relate the integrated MBH mass density to the integrated emissivity 



of the AGN population, via the integral of the luminosity function of quasars. If 
the average efficiency of converting accreted mass into luminosity is e = Lj Mc 2 , 
then the MBH will increase its mass by Mbh = (1 — accounting for the 

fraction of the incoming mass that is radiated away. Applying this argument 
to the whole MBH population, the MBH mass density can be related to the 
integral of the LF of quasar, *&(L,z), with the radiative efficiency being a free 
parameter. 

1.2. Relativist ic jets 

The so-called "spin paradigm" asserts that p owerful relativistic jet s are pro- 
duced in AGNs with fast rotating black holes ( Blandford et al. 19901 ). implying 



that MBHs rotate slowly in radio-q uiet quasars, which represent the majority 
of quasars ( Wilson fe Colbert 19951 ). However, if we expect the same mathe- 



matical and physical properties to describe both the black holes of a few solar 
masses and MBHs, such conjecture, at least in its basic inter pretation, is at 
odds with studies of radi o-emission from X-ray binaries (e.g.. lUlvestad fc Hoi 
2001 : Kording et al. 2006). These works showed that the production of jets is 



intermittent, although the spin of stellar mass black holes is not expected to 
vary over the same timescales. Why and when jet form, and what is the role of 
spin (if any) in jet production has not been explained yet. 

1.3. Gravitational recoil 

MBH spins affect the frequency of MBHs in galaxies, via the "gravitational re- 
coil" mechanism. When the members of a black hole binary coalesce, the center 
of mass of the coalescing system recoils due to the non-zero net linear momen- 
tum carried away by gravitational waves in the coalescence. If this recoil were 
sufficiently violent, the merged hole would breaks loose from the host and leave 
an empty nest. Recent breakthroughs in numerical relativity have allowed reli- 
able computations of black hole mergers and recoil velocities, taking the effects 
of spin into account. Non-spinning MBHs, or binaries where MBH spins are 
aligned with the orbital angular momentum are expected to recoil with veloci- 
ties below 200 kms . The recoil is much larger, up t o thousands kms" 1 , for 



MBH s with large spins in non-aligned co nfigurations (jCampanelli et al.l 120071 : 
Gonzalez et~aill2007l : iHerrmann et alj|2007l ). 
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2. Cosmic evolution of MBH spins 

MBH spins determine directly the radiative efficiency of quasars. On the other 
hand, accretion determines MBH spins. A hole that is initially non-rotating 
gets spun up to a maximally-rotating state (a = 1) after increasing its mass 
by a factor y6 — 2.4. A maximally-rotating hole is spun down by retrograde 
accretion to a = after growing by a factor ^/3/2 ~ 1.22. A 180° flip of the spin 
of an extreme-Kerr hole will occur after tripling its mass. Spin-up is therefore a 
natural consequence of prolonged disk-mode accretion: any hole that increases 
substantially its mass by capturing material with constant angular momentum 
axis would ends up spinning rapidly ("coherent accretion"). If the lifetime of 
quasars is long enough that angular momentum coupling between black holes 
and accretion discs through the Bardeen-Petterson effect effectively forces the 
inner most region of a c cretio n discs to align, then quasar MBHs should have large 
spins Volonteri et al. ( 20051 ). Spin-down occurs when counter-rotating material 



is accreted, if the angular momentum of the accretion disk is strongly misaligned 
with respect to the direction of the MBH spin. If accre tion proceeds via sma ll 
(and short) uncorrelated episodes ("chaotic accretion", King fc Pringld 12006). 



where accretion of co-rotating and counter-rotating material is equally probable 
then spins tend to be low. This is because counter-rotating material spins MBHs 
down more efficiently than co-rotating material spins them up (as the ISCO for 
a retrograde orbit is at larger radii than for a prograde orbit, the transfer of 
angular momentum is more efficient in the former case.) 



MBH-MBH mergers also influence the spin evolution. iBerti Sz Volonteri 



(2008) consider how the dynamics of BH mergers influences the final spin. Ex- 



cept in the case of aligned mergers, a sequence of BH mergers can lead to large 
spins > 0.9 only if MBHs start already with large spins and they do not expe- 
rience many major mergers. Therefore, the common assumption that mergers 
between MBHs of similar mass always lead to large spins needs to be revised. 



3. Spin evolution in gas-rich merger remnants 

Gas-rich mergers between galaxies of comparable mass (i.e. major mergers), are 
advocated to drive quasar activity, hence MBH growth. The evolution of MBH 
spins has to be addressed by simulating the typical environmental conditions ex- 
pected around quasars, where the properties of the accretion flow are established. 
When accretion is triggered by galaxy mergers, as expected for quasars, the ma- 
terial that will feed MBHs is expected to assemble into a circumnuclear disk. 
These disks may be the end-product of gas-dynamical, gravitational torques 
excited during t he merger, when large amounts of gas are driven into the core 
of the remnant ( Maver et al.ll2007l ). This last phase of the merger has been sug- 



gested to be associated to quasar activity. These are therefore the most relevant 
circumstances to be explored. During their inspiral MBHs are surrounded by 
a dense cocoon of gas t hat drives their dynami cal decay and provides fuel for 
the feeding of the holes dDotti et alJl200d . I2007D . Since matter carries angular 



momentum also the spin vector can change during the accretion process. The 
details of the dynamics may have a profound influence on the mass and spin 
evolution of the two MBHs. 
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3.1. Simulation set-up 

Dotti et al.l (|2010l ) follow the dynamics of MBH pairs in nuclear discs using 



numerical simulations run with the N-Body/SPH code GADGE T (Springel, 



Yosh ida & White 2001), upgraded to include the accretion physics (jDotti et al 



2009). In our models, two MBHs are placed in the plane of a massive circumnu- 



clear gaseous disc, embedded in a larger stellar spheroid. The disk is modeled 
with 2 xlO 6 particles, has a total mass 10 8 Mq, and follows a Mestel surface 
density profile. The spatial resolution of the hydrodynamical force is ~ 0.1 pc. 
With this spatial resolution the Bondi radius of the holes is resolved. The disk 
is rotationally supported and Toomre stable. SPH disk particles evolve with a 
purely adiabatic equation of state, described by either a polytropic in dex 7 = 7/5 
that mimics a star-forming region ("cold" disk. lSpaans fc SilTO OO). or 7 = 5/3 



that includes extra heat, e.g, AGN feedback ("hot" disk; see, e.g. iMaver et alj 
20071 ). 



The spheroidal component (bulge) is modeled with 10 5 collisionless parti- 
cles, initially distributed as a Plummer sphere with a total mass ~ 7x the disc 
mass. The mass of the bu lge within 100 pc is five times the mass of the disc 
(jDownes k Solomon! [19981 ) . The two MBHs (M 1 and M 2 ) are equal m mass 
(M BH = 4 x 10 6 M ). Mi, called primary for reference, is placed at rest at 
the centre of the circumnuclear disc. M 2 , termed secondary, is moving on an 
initially eccentric (eo — 0.7) counterrotating (retrograde MBH) or corotating 
(prograde MBH) orbit with respect to the circumnuclear disc. Gas particles are 
accreted onto a MBH if they lie within its Bondi radius, and if the total mass 
accreted onto a MBH in a timestep is lower than Mgdd- We find accretion rates 
close to Eddington for the central MBH (Mi) is in good agreement with the 
suggestion that quasars shine during the last phases of galaxy mergers. The 
secondary MBH in the simulation (M2) accretes instead at sub-Eddington rates 
during most of its orbital decay. 

Each gas particle accreted by the MBH carries with it angular momentum. 
From the properties of the accreted particles we can compute, as a function of 
time, the mass accretion rate and the versor l e dge> that defines the direction of 
the total angular momentum of the accreted particles. This information can be 
gathered only by performing very high resolution simulations. In our runs, the 
spatial resolution of the hydrodynamical force in the highest density regions is 
~ 0.1 pc, well below the Bondi radius. 



3.2. Modelling the evolution of spin vectors 

We use the MBH accretion histories obtained from our SPH simulations to follow 
the evolution of each MBH spin vector, Jbh = (aGM^ K /c) Jbh, where < a < 1 
is the adimensional spin parameter and Jbh is the spin versor. The scheme we 
adopt to study the spin evolution is based on the model recently developed by 
Perego et al. (2009). Here we summarize the algorithm used. 

We assume that during any accretion event recorded in our SPH simulations, 
the inflowing gas forms a geometrically thin/optically thick a-disc (Shakura & 
Sunyaev 1973) on milli-parsec scales (not resolved in the simulation), and that 
the outer disc orientation is defined by the unit vector l e d ge . The evolution of the 
a-disc is related to the radial viscosity v\ and the vertical viscosity v 2 : v\ is the 
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standard radial shear viscosity while v<i is the vertical shear viscosity associated 
to the diffusion of vertical warps through the disc. The two viscosities can be 
described in terms of two different dimensionless viscosity parameters, a\ and 
«2) through the relations = ax^Hcg, where H is the disc vertical scale 
height and c s is the sound speed of the gas in the accretion disc. We further 
assume a 2 = /2/(2ai), with a\ = 0.1 and f 2 = 0.6 (Lodato & Pringle 2007). We 
assume power law profiles for the two viscosities, v\ 2 oc as in the Shakura 

h Sunyaev solution. 

As shown by Bardeen Sz Petterson (1975), if the orbital angular momentum 
of the disc around the MBH is misaligned with respect to the MBH spin, the 
coupled action of viscosity and relativistic Lense-Thirring precession warps the 
disc in its innermost region forcing the fluid to rotate in the equatorial plane 
of the spinning MBH. The timescale of propagation of the warp is short com- 
parted with the viscous/accretion timescale so that the deformed disc reaches 
an equilibrium profile that can be computed by solving for the equation 



1 d , \ Id/ „„,dn A 1 9 /l „ r 91 \ 2G J BH x L . , 
— — (RLv R ) = — — i/iSi? 3 — 1 + — — -u 2 RL— +—-^— — (2) 
RdR y K; RdR\ dR ) RdR\2 dRJ c 2 i? 3 v; 

where vr is the radial drift velocity, E is the surface density, and 0, is the 
Keplerian velocity of the gas in the disc. L is the local angular momentum 
surface density of the disc, defined by its modulus L and the versor 1 that define 
its direction. 

The boundary conditions to eq. [2] are the direction of L at the outer edge 
ledge; the mass accretion rate (that fixes the magnitude of E), and the values of 
mass and spin of each MBH. All these values but the MBH spins are directly 
obtained from the SPH runs. In particular, the direction of the unit vector l e dge 
is computed considering those SPH particles nearing the MBH gravitational 
sphere of influence that are accreted according to the criteria outlined above. 

The MBH spin changes also in response to its gravito-magnetic interaction 
with the disc on a timescale longer than the time scale of warp propagation 
(Perego et al. 2009). This interaction tends to reduce the degree of misalignment 
between the disc and the spin decreasing with time the angle between Jbh and 
ledge- The MBH spin evolution is followed by solving the following equation: 

^ = MA(R lso )l(R lso ) + ^ [ ^^dR. (3) 
at c A ./disc R z 

The first term in eq. [3] accounts for the angular momentum deposited onto 
the MBH by the accreted particles at the innermost stable orbit (ISO), where 
A(-Riso) denotes the specific angular momentum at -Riso an d l(-Riso) the unit 
vector parallel to Jbh 5 describing the warped disc according to the Bardeen- 
Petterson effect. The second term instead accounts for the gravo- magnetic in- 
teraction of the MBH spin with the warped disc. It modifies only the MBH 
spin direction (and not its modulus), conserving the total angular momentum 
of the composite (MBH+disc) system (King et al. 2005). We applied iteratively 
eq. 1 and 2 using inputs from the SPH simulation that give the values of the 
mass accretion rate, the MBH mass and the direction of l e dge- The algorithm 
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Figure 1. Upper panels: time evolution of the relative angle between Mi 
(secondary MBH) spin and the orbital angular momentum of the MBH pair. 
Left (right) panel refers to a "cold" ( "hot" ) disk. The initial angle is arbitrar- 
ily set to 2.5 radians (close to anti-aligned), and the initial spin parameter 
magnitudes varies between 0.2 (lighter colours) to 1 (darker colours). Lower 
panels: same as upper panels for M\ (primary MBH). 
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Figure 2. Left: spin evolution for MBHs embedded in a "cold" circumnuclear 
disk. In a "cold" disk, where turbulence and pressure are relatively unimportant, 
most of the accreted material has angular momentum directed along the same axis, 
leading to mostly "coherent" accretion. Right: same quantities, but for a MBH 
evolving in a more chaotic "hot" circumnuclear disk. Here the direction of the 
angular momentum of the accreted gas is much more variable, giving rise to a 
"chaotic" flow, and the magnitude of the MBH spins ends up being lower, at a 
given time along the simulation. Bottom: primary, central, MBH. Top: secondary 
MBH. In each panel we show 100 realizations starting from a flat distribution in 
spin magnitudes and initial orientations. 
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returns, as output, the spin vector, that is, its magnitude and direction. At each 
timestep our code therefore provides the angle between the spin vector of each 
MBH and the angular momentum vector of their relative orbit. 

3.3. Spins: alignment and magnitudes 

Figure [1] shows the time evolution of the relative angle 9 between the spin of each 
MBH and the orbital angular momentum of the MBH pair. The initial relative 
angle has been arbitrarily set to 2.5 radians (143°), while a has initially five 
different values (0.2, 0.4, 0.6, 0.8, and 1). In all cases we find that MBH spins 
loose memory of their initial orientation: accretion torques suffice to align the 
spins with the angular momentum of their orbit on a short timescale (< 1 — 2 
Myr). A residual off-set in the spin direction relative to the orbital angular 
momentum remains, at the level of ^ 10° for the case of a cold disc, and ^ 30° 
for a warmer disc. Alignment in a cooler disc is more effective due to the 
higher coherence of the accretion flow near each MBH that reflects the large- 
scale coherence of the disc's rotation. If the MBHs coalesce preserving the spin 
directions set after formation of a Keplerian binary, the MBH resulting from 
their coalescence receives a relative small kick, < 100 km s -1 . 

Fig.[2]shows the time evolution of the spin magnitudes (Dotti et al. in prepa- 
ration). We show here 100 realizations starting from a flat distribution in spin 
magnitudes and initial orientations. We find that initially maximally rotating 
are spun down, and initially slowly rotating holes are spun up, leading to inter- 
mediate equilibrium values (~ 0.6 — 0.8). The spin of the secondary MBH tends 
to be slightly lower, because of a higher incidence of counterrotating events. It is 
suggestive that measured spins cluster around these values. The best measure- 
ment to date is the Seyfert galaxy MCG-6-30-15 . In this Seyfert, the iron line is 
so broad as to rigoro usly require a = 0.9 ±0.1 (|Brenneman & Reynolds! 120061 ) . 



Miniutti et aD (12009T ) reported a spin parameter of a = 0.6 ± 0.2 in the narrow- 



line Seyfert-1 AGN SWIFT J2127.4+5654. The X -ray spectrum of Fa irall 9 also 



results in a spin a = 0.7 ± 0.2 at 99% confidence (jSchmoll et al.ll2009l ) 



4. MBH spins and galaxy morphology 

If the events powering quasars coincide with the formation of elliptical galaxies 
(di Matteo et al. 2005), we might expect that the MBH hosted by an elliptical 
galaxy had, as last major accretion episode, a large increase in its mass. During 
this episode the spin increased significantly as well, possibly up to very high 
values, a ~ 0.6 — 0.8, as suggested in the previous section. 

Black holes in spiral galaxies, on the other hand, probably had their last 
major merger (i.e., last major accretion episode), if any, at high redshift, so 
that enough time elapsed for the galaxy disc to reform. Moreover, several 
observations suggest that single accretion events last ~ 10 5 years in Seyfert 
galaxies, while the total activity lifetime (b ased on the fraction of disc galax - 
ies that are Seyfert) is 10 8 - 10 9 years fe.g-. lKharb et al.ll2006l : iHo et al.lll997h . 



This suggests that accretion events are very small and very 'compact'. Smaller 
MBHs, powering low luminosity AGN, likely grow by accreting smaller packets 

1 , 



of material, such as tidally disrupted stars (for MBHs with mass < 2 x 10 6 M^ 
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Milosavlievic et all [20 06). or possibly molecular clouds ( Hopkins Sz Hernquist 
20061). 



Compact self-gravitating cores of molecular clouds (MC) can occasionally 
reach subparsec regions. Although the rate of such events is uncertain, we 
can adopt the estimates of Kharb et al. (2006), and assume that about 10 4 
of such events happen within the total activity lifetime of a Seyfert. We can 
further assume a lognormal distribution (peaked at log(MMc/ M Q ) = 4, with a 
dispersion of 0.75) for the m ass function of MC close to galaxy centers (based 
on the Milky Way case, e.g., iPerets et al.1 120071 ). The result is, on the whole, 
simil ar to that produced by minor mergers of black h oles (Hughes &: B landford 



20031 1. that is a spin down in a random walk fashion (jVolonteri et al.l [2007) 



In a gas-poor elliptical galaxy, however, substantial populations of molecular 
clouds are lacking (e.g., Sage et al. 2007), eliminating this channel of MBH 
feeding. Main sequence stars, however, linger in galaxy centers. Tida l disruption 
of stars is a feeding mechanism that has been proposed long ago ( Hills! 19751 ; 



Reed 19881 ) . One expects discs formed by stellar debris to form with a random 



orientation. Stellar disruptions would therefore contribute to the spin-down 
of MBHs. In an isothermal cusp, assuming that MBH masses scale with the 
velocity dispersion, a, of the galaxy (we adopt here the Tremaine et al. 2002 
scaling), we can derive the relative mass increase for a MBH in 1 billion years: 

m * o, 7 (j^r /8 . (4) 



M BH ' V1O 6 M 0/ 

The maximal level of spin down would occur assuming that all the tidal dis- 
ruption events form counterrotating discs, leading to retrograde accretion. Eq. 
H] shows that a small (say 10 5 M Q ) MBH starting at a = 0.998 would be spun 
down completely, as its mass increase is larger than \/3/2 (cfr. section 2). On 
the other hand the spin of a larger (say 1O 7 M0) MBH would not be changed 
drastically. This feeding channel is likely efficient in early type discs which typ- 
ically host faint bulges characterized by steep cusps, as exemplified above. The 
situation is different for giant ellipticals: the central density profile displays a 
shal low core, and tidal dis ruption of stars is unlikely. 

Volonteri et a l. (2007) therefore suggest the spin of MBHs in giant elliptical 



galaxies is likely dominated by massive accretion events which follow gas-rich 
galaxy mergers (see Dotti et al. in this volume for a discussion of spin evolution 
in ellipticals formed in gas-poor mergers). Both tidal disruption of stars, and 
accretion of gaseous clouds is unlikely in shallow, stellar dominated galaxy cores. 
In a galaxy displaying instead power-law (cuspy) brightness profiles, the rate of 
stellar tidal disruptions is much higher and random small mass accretion events 
contribute to spin MBHs down. This result is in agreement with Sikora et al. 
(2007) who found that disc galaxies tend to be weaker radio sources with respect 
to elliptical hosts. 
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